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■ Abstract 

We study the stability of layered structures in a variational model for diblock copolymer- 
(— I ' homopolymer blends with respect to perturbations of their interfaces. The main step consists 

^ ^, of calculating the first and second derivatives of a sharp-interface Ohta-Kawasaki energy for 

I ■ straight mono- and bilayers and determining the latter's sign. By developing the interface 

i^H ' perturbations in a Fourier series we fully characterise the stability of the structures in terms 

, of the energy parameters. Both for the monolayer and for the bilayer there exist parameter 

regions where these structures are unstable. For strong repulsive interaction between the 
monomer types in the diblock copolymer the bilayer is always stable with respect to interface 
perturbations, irrespective of the domain size. The monolayer is only stable for small domain 
size. 

^ I In the course of our computations we also give a Green's function for the Laplacian on a 

. two-dimensional periodic strip. 

0\ . 
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' 1 Introduction 



1.1 Localised and partially localised patterns 

Localised patterns are observed in a wide variety of systems, including experimental systems 
such as the Belusov-Zabotinsky reaction [5S], nonUnear optics |48[ HS] . vertically shaken granular 
media [5^ [5D] , and Bose-Einstein condensates [J^ i a-nd also in idealised systems such as the Swift- 
Hohenberg equation [Til HH [45l [49] or networks of reacting cells [29]. More recently objects have 
been observed that are only partially localised; structures in two dimensions, for instance, that are 
'thin' in one spatial direction and 'long' in the other. Such partially localised patterns have been 
observed in Nonlinear Schrodinger equations [131 |S1|T1|51|3] , Gierer-Meinhardt-type systems [Tl] , 
and even in scalar nonlinear elliptic equations [261 1271 125] . In addition, the membrane that 
surrounds each living cell, for instance, is such a structure [23l [7l [34]. 

In this paper we study an example of energy- driven partial localisation, arising in the study of 
mixtures of diblock copolymers with homopolymers. Such mixtures feature two opposing forces: a 
repelling force between different monomer types favours separation into homogeneous phases, while 
covalcnt bonds between some of the repelling monomers impose an upper limit on the separation 
length. As a result a wide variety of patterns are observed (both in physical and in numerical 
experiments) , ranging from spheres [22l [33l [53l [57] , cylinders [21] , dumbbells [32] , helices [19] , 
'labyrinths' and 'sponges' [H [201 [22, 'ball-of-thread' [11], layered structures [ ^ [^[5 ^ [57]. and 
many more. 
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Our focus is on layered patterns, consisting of two or more parallel layers of roughly uniform 
thickness. In each layer the composition is dominated by one of the polymer types, and in the 
separation into layers one can recognise a phase separation phenomenon triggered by the repelling 
forces between polymer types. In addition to their interest as particular patterns in copolymer- 
homopolymer blends, such layered structures are examples of energy-driven partial localisation. 

The main goal of this article is to understand the (in)stability of such layered structures in 
this simple model of copolymer- homopolymer blends. 



1.2 Diblock copolymers and blends 

Diblock copolymers are linear polymer molecules that consist of two parts (blocks) called the U- 
part and the V-part in this paper, with corresponding volume fractions given by the functions u 
and V. Each part contains monomers of a single type only, U or V. As described above, the 
interaction between the two types of monomers is the net result of two opposing influences. On 
the one hand the U- and V-parts repel each other, leading to a tendency of the U-and V-phases to 
separate; on the other hand the U-and V-parts are chemically bonded together in a single diblock 
copolymer molecule, forcing both parts to remain close to each other. As a result of these two 
types of interaction, the separation between the U-and V-phases is restricted to length scales of 
the order of the molecule size. 

We consider systems that contain, in addition to the diblock copolymers, some species of 
homopolymer, that we call the 0-phase. A homopolymer is made up of a single type of monomers, 
here named 0. The system therefore contains three phases, and because of an assumption of 
incompressibility we can use the functions u and v to describe the distributions of the three 
phases. 

In [8] the following energy is derived: 

I Co \V{u + v)\+Cu \Vu\+Cy \Vv\ + ||u-i'||^-i(St) if(",w)G^, 
■r(u,v) - Js^ Js^ Js^ 

oo otherwise, 

where the coefficients q are nonnegative (and not all equal to zero), Sl is a periodic strip x R 
(where is the one-dimensional torus of length L) , and the set of admissible functions is given 
by 

IC -.^ Uu,v) e {BY {SL)f ■ u{x),v{x) e {0,1} a,.c., uv = a..c., and / u = / v 

Since unconstrained minimisation will lead to the trivial structure u = v = 0, the natural problem 
to look at here is minimisation under constrained mass, i.e. with the constraint J^^ u = J^^ v ~ M 
for some M > 0. 

Under the extra restriction u -\- v = 1 — no 0-phase — the functional J- la a, well-known sharp- 
interface model for diblock copolymer melts |351 [9]. The sharp-interface character of this model, 
known in the physics literature as the strong-segregation limit, is recognizable in the fact that the 
variables u and v are characteristic functions, implying that at each point in space only one phase 
is present. The underlying diffuse-interface model is weU studied [3T | [T6 l fT5 l [50 l [36 l [9l[37 l l4g i [4T | 
142 1 15 H l43] because of the interesting pattern formation phenomena it exhibits. 

The first three terms of T can be recognised as the sharp- interface manifestation of the repelling 
forces between the U-, V-, and 0-monomers. The last term, the i/~^-norm, is a remainder of the 
chemical bond between the U- and V-parts and penalises large-scale separation of the U- and 
V-phases. 

The functional J- resembles the energy functional used to model triblock copolymers, i.e. block 
copolymers consisting of three chemically bonded parts, [38, 39J. The interface penalisation part 
is present in that functional as well, and the long range interaction term includes interaction 
between the third phase (the phase corresponding to the third part of the triblock copolymers) 
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and the other two phases in addition to the interaction between the first two phases present in 
the functional T above. 

For a more extensive review of the modehing of dibloek eopolymers and dibloek copolymer- 
homopolymer blends and its study in mathematics, we refer to [541 Chapter 2]. 

1.3 From one-dimensional to two-dimensional structures 

A layered structure with perfectly straight layers can be described by functions u and v of one 
spatial variable. In a companion paper [17j (see also [8]) we study this one-dimensional case and 
give a full characterisation of global minimisers. 

One of the results in that paper is that, for generic parameter values, every constrained-mass 
global minimiser on M is a concatenation of equal-width monolayers. A monolayer is shown in 
Figure [TJ a structure, described by a pair of functions (it, v), in which the supports of u and v are 
adjacent intervals of equal length — or, in the higher-dimensional context, adjacent layers of equal 
width (see Figures [Tal and [Tc)) . 

For small constrained mass, the global minimiser in one dimension is a monolayer. For slightly 
larger constrained mass, the global minimiser switches to a bilayer, a pair of monolayers joined 
back to back (Figures llbl and [Td)) . As the constrained mass further increases the global minimiser 
switches to structures of increasing numbers of monolayers (see [H])- 

In the present paper we are interested in the stability properties under of a particular subset 
of two-dimensional mono- and bilayer structures: 

• For both mono- and bilayers we assume that the layer thickness is such that the energy-to- 
mass ratio J- / J u is minimal among all such layers; 

• For monolayers we assume that c„ = c„, i.e. that the interface penalisation is the same for 
U-0 and V-0 interfaces. 

Both restrictions arise from our interest in thin, partially localised structures in R^, as is explained 
in detail in Appendices El and [BJ 

The optimal widths Sm and St for which the energy-to-mass ratio is minimal for the mono- 
and bilayer respectively are indicated in Figure [T] and defined in pT]) and ([7]). 

1.4 Stability of mono- and bilayers in two dimensions 

The aim of this paper is to investigate the stability of these mono- and bilayers in two dimensions. 
Since the functions u and v are forced to be characteristic functions of sets, the only admissible 
perturbations are changes in the supports of these functions. In this paper we only consider local 
stability with respect to perturbation of the position of the interfaces; other perturbations, such as 
those that change the topology of the structure are disregarded (see the discussion in Section [6]). 

Specifically, we consider perturbations of the interfaces that are periodic with period L along 
the length of the layer, and therefore we assume a domain that is periodic in one direction (xi) 
and unbounded in the other (see Figure [T|) . Because of this periodicity each perturbation of an 
interface is given by a periodic function p : Tl ^ M.^ (for the monolayer) or p : ^ (for the 
bilayer), where each component is the lateral displacement of one of the interfaces. By expanding 
the perturbations in Fourier modes, and using the usual vanishing of cross terms of different 
frequency, thcpositivity of the second derivative of the energy reduces to the positivity on each 
Fourier modcH 




Fourier modes have a natural scale invariance: the k Fourier mode on the interval of length L 
is equivalent to the Fourier mode on an interval of length L/k. This allows us to establish 
the stability with respect to the first Fourier mode as a function of L, rescale for the stability 
properties of the k*"^ mode, and aggregate the results. 

^ First and second derivatives of similar functionals have been calculated by Muratov and Choksi & Sternberg 1301 
1101 . Our calculations differ in the number of phases (three instead of two) and in the early adoption of a Fourier 
framework. 
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(a) A one-dimcnsional monolayer 
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(c) A straight monolayer on the periodic 
strip Sl 
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(b) A one-dimensional bilayer 
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(d) A straight VUV bilayer on the periodic 
strip Si 



Figure 1: Mono- and bilaycrs on a strip as trivial extensions of one-dimensional structures. We 
assume the O-phase to fill up the rest of the domain where there is no U- or V-phase. We do not 
indicate this in our pictures. 



Using this approach we show in Section |3] that the monolayer of optimal width is linearly 
stable with respect to mode-1 perturbations iff 



> fliL/5rn), 



2c„ + Co 

where /i is an explicit function given in (|3ip . By combining all Fourier modes we find 
Theorem 1.1. Assume Cu = Cy. The monolayer of optimal width is linearly stable iff 



2cu + Co 



> f{L/5m,) 



where 



/(£) :=sup/i(^/A:). 

k>l 



The graphs of the functions x i-^ fi{x/k) for different values of k are shown in Figure !^ 

For a bilayer of optimal width 5h a similar result holds: 
Theorem 1.2. The VUV bilayer of optimal width is linearly stable iff 



Co + C„ -I- 2Cy 



> 9{L/h) 



where 



and gi is given by 



g[i) := sup5iWA:) 



(1) 



(2) 
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/i(L/(M™)) 



9i{Ll{k5,)) 
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(a) The monolayer; plotted are the curves L/(5, 
/l(L/(fc5™)) for fe = 1...20 



(b) The bilayer; plotted are the curves -L/<5(, 
gi{L/{k5i)) for A: = 1. ..20 



Figure 2: The graphs of the funetions x fi{x/k) and x i— > gi(x/k) {k ~ 1, . . . , 20) portray the 
curves in parameter space that separate the parts where the first twenty Fourier modes of the 
second variation for the monolayer (Figure [2a|) and bilayer (Figure [2b)) are positive and negative. 
If c„/(2c„ + cq) < fi{L/{kSm)) the fcth Fourier mode is negative for the monolayer, if the reverse 
inequality holds the mode is positive. Similarly for the bilayer the fcth Fourier mode is negative 
if {cu + Cv)/{co + Cu + 2c„) < (?i(L/(fc(S(,)). The leftmost curve in each figure corresponds to the 
first order Fourier mode, the order increases towards the right. Note that the positivity of the 
parameters c„ and cq implies that c„/(2c„ + co) < 5 as indicated in Figure [2al bv the dashed line. 

In Figure [2bl the graphs of the functions x 1-^ gi{x/k) are shown for different values of fc. 

From Figures [^a] and I2bl one might think that curves belonging to higher orders remain below 
curves of lower orders. The blow-ups in Figure [3] however show that this is not the case. However, 
it is true that only the first Fourier mode is of importance for determining the stability of the 
monolayer. This can be recognised by noting that the left-hand side of ^ cannot reach values 
larger than 1/2, and Figures [2al and [3al show that the non-monotonicity for the monolayer plays 
a role only for values above 1/2. 



Figure 3: A blow-up of the graphs in Figure [H Curves corresponding to different Fourier modes 
clearly cross. Here ^ = 2(eo+^°+c ) ^'^'^ ~ eo+c^+2c ' ^^^"^ ([^0)) and (f5(Il) . 

Figure |4] summarises the stability properties of both the mono- and the bilayer. In Figure l4al 
the vertical axis is restricted to the interval [0, 1/2] to reflect the value set of the left-hand side 
of ll]). This implies that monolayers can only be stable if L is sufficiently small, and even then 
only for a subset of the coefficients cq, c„, and c^,; for sufficiently large L the monolayer is unstable 
for all choices of interface penalisation. 

For the bilayer the situation is different: here the condition ([2]) allows for both stability and 
instability at all values of L. The function g is bounded from above (away from 1), implying that 




(a) Blow-up of Figure I2al 



(b) Blow-up of Figure I2bl 
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a threshold a exists such that 




V 



0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 

(a) Monolayer (b) Bilayer 



Figure 4: The sign of the second derivative operator for the mono- and bilayer of optimal width. 
+/— indicates indeterminate sign, due to the negativity of one or more eigenvalues. Along the 
horizontal axes are plotted z/ = g-27r5„/L _ ^-2Tr5b/L^ rpj_^g vertical axes show /i = c„/(2c„ + 

Co) and ( = (c„ + Cy)/{co + c„ + 2c„). These figures are based on a calculation involving Fourier 
modes up to and including order 100. 



1.5 Directions of instability 

For the functional J- one may imagine a number of different evolution problems, such as gradient 
flows based on the L^, H^^, or Wasserstein metrics. Under such an evolution the straight mono- 
and bilayer structures are stationary. If they are unstable, the evolution will amplify small devi- 
ations and move away from the straight configurations. While the perturbations are still small, 
the main contribution of the evolution will be in the directions of the eigenvectors of the second 
variatiorll belonging to the (most) negative eigenvalues. 

For the monolayer there is, for each Fourier mode, one eigenvalue that can become negative 
(for the first Fourier mode: in Lemma r4.12( other modes follow by rescaling as above) and there 
are two which are always positive. Each component of the corresponding eigenvectors is associated 
with the deformation of one of the interfaces in the layer. A cartoon of the (possibly) unstable 
deformation direction is given in Figure ISa) the two stable directions are shown in Figures I5bl 
and Eel 

For the bilayer two eigenvalues are always positive, and two eigenvalues may also become 
negative. For the first Fourier mode the dependence of the sign of the latter two on the parameters 
L /6b and ( is given in Figures [6a| and I6bl We recognise in the second figure the first order curve 
(fc = 1) from Figure [2bl a similar curve for the first figure would always stay below the curve from 
the latter one, which is why its influence is not recognisable in Figure [2bl 

The (possibly) unstable deformation directions are shown in Figures iTal and Ffbl corresponding 
to the eigenvalues in Figures [Hal and I6bl the stable ones in Figures [7c] and I7dl 

These results all show that depending on the parameters in the model the monolayer and 
bilayer structures can be unstable. This mirrors closely the results in [30], [37], and [40], where 
it is shown that in the pure diblock case 'wriggled' lamellar structures bifurcate off the straight 

^For each Fourier mode the bihnear form that is the second variation can be identified with a bilinear form 
on M.^ (monolayer) or (bilayer) whose eigenvalues and eigenvectors can be studied. Details can be found in 
Sections |42HO] 
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(a) The (possibly) un- (b) One of the stable di- 

stable deformation direc- rections 

tion 




(c) The other stable di- 
rection 



Figure 5: One (possibly) unstable and two stable first order Fourier modes of deformation for the 
monolayer; see Section |4l in particular Remark [4. 171 



c c 




(a) The sign in parameter space of the (b) The sign in parameter space of the 

eigenvalue corresponding to the eigenvalue eigenvalue corresponding to the eigenvalue 

G-i- of the reduced matrix Bi in the proof G- of the reduced matrix Bi in the proof 

of Lemma 14.51 of Lcmma l4.5l 



Figure 6: The black patches in parameter space indicate where two of the eigenvalues of the first 
Fourier order second variation operator for the bilayer become negative. 

lamellar pattern if the spacing between the lamellae becomes too large. In Section [6?3l we discuss 
the relation with these results in more detail. 



1.6 Structure of this paper 

We start in Section [5] by defining the functional under consideration and clarifying some of the 
notation that is used throughout the paper. In Section [3] we prove via a calculation of the first 
variation of that the monolayer and bilayer are both stationary points of J- with respect to mass 
preserving perturbations of the interfaces. We then proceed to compute the second variation for 
both these structures. Since this calculation for the monolayer is similar to that for the bilayer, 
we only give the details in the latter case and even there we defer most of the computational 
details to Appendix [Cl Section |4] is dedicated to computing the sign of the second variations for 
the monolayer and bilayer in order to determine the parameter regions of stability and instability. 
Much of the work in the proofs is again of a calculational nature, some of which we have also 
moved to the back of the paper in Appendix [P] Finally Section [5] gives a Green's function of — A 
on the periodic strip. This Green's function is heavily used in this paper and since the authors 
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(a) One of the (possibly) 
unstable deformation di- 
rections 











u 







(c) One of the stable di- 
rections 



(b) The other (possibly) 
unstable deformation di- 
rection 









u 





(d) The other stable di- 
rection 



Figure 7: Two (possibly) unstable and two stable first-order Fourier modes of deformation for 
the bilayer. The (possibly) unstable deformation in Figure [7a| corresponds to the eigenvalue in 
Figure [6a| and the deformation in Figure [7b] to the eigenvalue in Figure [6bl For details see the 
discussion in Section IH in particular Remark 14.101 



could not trace a previous appearance of it in the literature a section on its validity closes this 
paper. 



2 Definitions and conventions 
2.1 Problem setting 

The domain of definition is the strip Sl ■= x M, where is the one-dimensional torus of 
length L, i.e. the interval [0, L] with the endpoints identified. For functions on Sl the i?~^-norm 
is defined by convolution: 

Definition 2.1. For f e L°°{Sl) with compact support satisfying J^^ f — we define 

11/11^-1(5^):= / / f{xi,X2)G * f{xi,X2)dx2dxi, 
Jq Js. 

where G is the Green's function of the operator —A on Sl, i-e. it satisfies —AG = 6 in the sense 
of distributions (6 is the Dirac delta distribution). 

Note that 0/ -.^ G * f satisfies — A^/ = f on Sl- Also note that while the Green's function 
is only unique up to addition of an affine function of X2, this non- uniqueness is irrelevant for the 
definition above. 

We repeat the definition of J- and IC for convenience. 

Definition 2.2. Let cq, c„, and Cy be real numbers. Define 



T{u, v) — 




if {u,v) e IC, 
otherwise. 



where the admissible set is given by 

K. -.^ l^{u,v) e {BV{SL)f ■■ u{x),v{x) e {0,1}, uv = a.e., and J u = J w| . 



Wc will require that all Ci are non-negative and at least one of them is positive. 
Another, equivalent, form of the functional will be useful, in which the penalisation of the three 
types of interface U-0, V-0, and U-V, is given explicitly by surface tension coefficients dkV- 



Lemma 2.3. Let the surface tension coefficients be given by 




Non-negativity of the Ci is equivalent to the conditions 



< dki < dkj + dji 



for each k ^ I ^ j ^ k. 



(3) 



Then 



J^{u, v) — 



duon^-\Suo) + dMn^-\S^o) + duvH^-^iSuv) + \\u - v\ 



if {u,v) e K,, 



oo 



otherwise. 



where Ski is the interface between the phases k and I: 



Suo = d* supp u\d* supp V, 
Svo = d* supp v\d* supp u, 
Suv = d* supp uC\d* supp V, 



and d* is the essential boundary of a set. 

The essential boundary of a set consists of all points in the set that have a density other than 
or 1 in the set; see e.g. [H Chapter 3.5]. 

Proof of Lemma \2.3[ The main step in recognising the equivalence of both forms of is noticing 
that, for characteristic functions of a set, such as u, v and u + v, the equality 



Note the different interpretations of the coefficients Ci and the surface tension coefficients 
dki- The latter have a direct physical interpretation (and can be related to material parameters, 
see j8]): they determine the mutual repulsion between the different constituents of the diblock 
copolymer-homopolymer blend. For example, the value of duv (as compared to the values of duo-, 
dyo and 1, the coefficient in front of the _ff~^-norm) determines the energy penalty associated with 
close proximity of U- and V-polymers. In particular, if one of these surface tension coefficients 
is zero, the corresponding polymers do not repel each other and many interfaces between their 
respective phases in the model can be expected. On the other hand the coefficients c^, when taken 
separately, do not convey complete information about the penalisation of the boundary of a phase. 
If for instance c„ = 0, but Cy 0, the part of the U-phase interface that borders on the V-phase 
still receives a penalty, because duv = Cv For this reason the use of surface tension coefficients 
makes more sense from a physical point of view. For the mathematics it is often easier to use the 
formulation in terms of Ci . 

If we consider the functional on three-dimensional 'physical space' we can also see from 
dimensional considerations that the name "surface tension coefficients" for the dki is justified. 
Since Ti.'^{Ski) measures surface area, if J- has the dimension of energy then the coefficients dki 
have the dimension of energy per unit area, which is also the dimension of surface tension. 

^The indices j, k, I take values in {u, v, 0} and the df^i are taken symmetric in their indices, i.e. dyu = d^v etc. 




(see [H Theorem 4.4], [1 Theorems 3.59, 3.61]). 



□ 
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Remark 2.4. The conditions ^ can be understood in several ways. If, for instance, d^v > 
duO + dyo, then the U-V type interface, which is penahsed with a weight of duv, is unstable, for 
the energy can be reduced by slightly separating the U and V regions and creating a thin zone 
of inbetween. A different way of seeing the necessity of ^ is by remarking that the equivalent 
requirement of non-negativity of the Ci is necessary for to be lower semicontinuous in e.g. the 
topology (see e.g. [HI Theorem 1.9]). Our assumption that at least one q is positive is equivalent 
to requiring at least two dki to be positive. 

2.2 Fourier transformation 

To clarify the notation we use, we will explicitly define the Fourier series we are using. For future 
reference we will also state some results we will need. 

Definition 2.5. Let f e L'^{Tl), then we will denote by f E _L'^(Z;C), the Fourier transform 
of /; 

V L Jo 

and by aj and bj , j G N, the Fourier coefficients of f with respect to the normalised basis of cosines 
and sines: 



ao ■- I f{x)dx, 



^ /(x)cos (^^^^ dx, 
bj := fix) sin (^^^^^ dx, 



For easy reference we give here the relations between f{j) and aj, bj: /(O) = ao and, for j > 1, 

/O') = 71 («^- - ^^^O' /(-•?) = 71 + *^^)' «^ = 71 (/(.?) + /(-.?)) and b, = ^ (/(j) - /(-.?)) ■ 
Furthermore we have 



ao h 2 
fix) ^ -Tf + \] j^^o-j^oa i27rxj/L) + J j^^bj sin {2nxj/L) 

.i=i .7=1 



oo 



where the convergence is in the topology. Finally Parseval's theorem takes the form 







fix)gix) dx = /(O)g(O) + 2Re ^ f{q)giq) 



q=l 

OO 



and as a consequence we have for pi,p2,P3 G L^(Tl) 

/ / Piix)p2ix-y)p3iy)dxdy = L^/^Ypiiq)p2iq)p3iq)- (5) 
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3 Geometrical derivatives of the energy 



In the following two sections we will take a look at the stability of two-dimensional periodic 
monolayer and bilayer configurations. First we need to determine under which conditions these 
structures are stationary points of the functional T . For the bilayer this will be done in Section [3. 11 
after which we compute the second variation for a bilayer in Section [3T2l We will give analogous 
results for the monolayer in Section [3731 In Section [4] we will use these results to derive the explicit 
stability criteria of Theorems 11.11 and 11.21 

Of the two possible bilayer structures — UVU and VUV — wc only discuss the VUV structure. 
The results for the UVU structure follow from exchanging the roles of u and v. 



3.1 Bilayer: admissible perturbations and stationarity 

The VUV bilayer of optimal width is a structure given by functions (mo,wo) G with 

where [T^ 

<56 := Y^^(co + c„ + 2c„), (7) 

and XA is the characteristic function of the set A. The set of admissible boundary perturbations 
of this structure is only restricted by regularity and the equal-mass constraint: 

Definition 3.1. The set of admissible perturbations is characterised by 

Vb ■= \p e {W^'\TL)t : 2 / {Pi+Pz)= I {P2 +P4)] . (8) 
For p Cz Vb and e > we define a perturbed structure (u^, v^), 

u (xi X2) = I ^ ^ {-6b - ep3{xi),db + epi{xi)) , 

^ 1 otherwise, 

^, ^ if X2 e {-2Sb ~ ep4,{xi),~Sb ~ ep3{xi)) U (5b + spi{xi),2Sb + ep2{xi)), 

^ ' [0 otherwise. 

We also introduce the subset of perturbations that conserve mass: 



V, 



M 



peVb-. {Pl+P3)^ {P2+Pa)^0 



(9) 



Note that since W^''^{Tl) is imbedded in L°°{Tl) by the Sobolev imbedding theorem, the pair 
belongs to JC for sufficiently small e. 
A picture of a bilayer of optimal width with perturbations p is shown in Figure [51 

2Sb + ep2{xi) 

Sb +epi{xi) 



■■ -5b - epsixi) 
25b - £P4ixi) 




Figure 8: The bilayer of optimal width with perturbations 
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Remark 3.2. We should stress the difference between the two mass constraints ^ and The 
constraint ([S]) is equivalent to the condition that and have the same mass. This property is 
a basic element of the model of block copolymers, via the set of admissible functions /C. 

The additional condition ^ expresses the requirement that J and J both equal the mass 
J uq of the unperturbed bilayer; perturbations without this property are meaningful in a situation 
where the joint mass of and may change. The bilayer of optimal width is a stationary point 
of the functional T under mass-preserving changes (see Lemma 13.41 below) : but as can be inferred 
from equation ()13p . the functional is not stationary under perturbations that do change the mass. 



Definition 3.3. We say that the bilayer of optimal width is stationary with respect to the admis- 
sible perturbations Vt (or V^^ ) if, for all p E Vb (or all p E V^^ ), 

d 



de 



= 0. 



Lemma 3.4. The VUV bilayer of optimal width is stationary with respect to all p G V^^ . 

Proof. Choksi and Sternberg calculate the first and second variations of a related functional [10] , 
and their method can be adapted without much difficulty to the functional J-. Here we give a 
self-contained proof. 

Since the interfaces of the bilayer are straight, the derivative of the interfacial terms with 
respect to the perturbation is zero for all p E Vb- 



d_ 




de 




d 




de 


= 0. 





I \Vu,\- 









e=0 

f (^^J 1 + ffVi' + \/l + f'Ps") dx + d,o + ^Va' + \/l + fM") dx 

(10) 



For the derivative of the H ^-norm, let i] E C(M) and compute 



— / r]{x2)ue{x) dx 



e=0 



de 



r]{x2) dx2 



e=0 



dxi 



(piixi) ri{Sb + epi{xi)) +P3{xi)r]i-Sb - ep3{xi))j dx^ 



e=0 



P3- 



(11) 



Similarly, 



d 
de 



£=0 







T]{x2)ve{x) dx ^ -rj{db) / pi + ri{2Sb) / p2-il{-5b) / + '?(-256) / Pa- (12) 











Let G be the Green's function of —A on Sl from Theorem 15.11 then 

6 = 

= 2 



d_ 

de 



■ dx 



6=0 



— VG * (Ue - v.] 
de 



f \VG*iu,-v,)\ 

JSl 

\ f 

VG * (uo - Wo) 
2 — / VG * (wo - vq) ■ \7G * (lie - v^ dx 

d£ J Sr. 

[G* (tio ~vo)]{u^ -Ve)dx 



dx 



6=0 



6=0 



de 



6 = 



12 



Setting ri{x2) [G* {uq — vq)]{xi, X2) (which is independent of xi, because uq — vq is independent 
of xi) we calculate by the Fourier series (j33p (or by remarking that this is a one-dimensional 
situation) that 

^7(2^2) = -^ J \x2-y\iuo-vo)i0,y) dy, 

from which it follows that ri{6b) = ri{—db) and i]{±2db) = 0. Therefore rj e C(R) and thus we 
obtain from ([Til) and ((12]) that 



AL^iSb) l{pi+P3)=0 



r=0 



(13) 



□ 



Remark 3.5. Note that in Lemma we nowhere use the specific definition of Sb in ([7]). As we 
explain in Appendix 1X1 the optimal width is relevant when considering energy per unit mass. If 
we define the mass functional as 



M{u, v) := I u 

'Sl 

for (u, v) G /C. then the following calculations show that the bilayer of optimal width is a stationary 
point oi T / Ai with respect to all perturbations in Vb (so not only the mass preserving ones) in 
the sense of Defintion 13.31 (with the functional T replaced by T j M). Thus, let now p G Vh- 

We first compute that "qiib) = \^\^ where 77 is such as chosen at the end of the proof of 
Lemma Then with the help of (fTS]) and the computations for the one-dimensional case in [T7j 
we find that 



d 
de 
d 
de 



e=0 



2<5fc' / (P1+P3), 



[pi +P3), 



E=0 

M{uo,vo) = 2LSb 



2L{ca + Cu + 2c^) + -LSI 



Now we conclude 

d_ 

de 



T/M{u„v,) 



e=0 



M{uo,vo) '^lM{uQ,vo)—J^ius,Vs) 



e=0 



de 



6=0/ 



\l 2 /" (pi + pg) i _ ^ ^ 2C„) 



0. 



We see that the optimal width condition ^ is not necessary for stationarity under T j M. with 
respect to the mass preserving perturbations in V^^ . but it is for stationarity with respect to 
perturbations in Vb \ V^' . 



3.2 Second variation for a bilayer 

We express the components pi of a given perturbation p e T'h as a Fourier series (see Section 12.2 



Pi{.x) 



7z 



^a,,cos(^-^j 



^ 6j J sir 



2-Kxi 
L 



(14) 
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The equal-mass condition in ([5]) translates into 

2 (ai,o + 03,0) = 02,0 + 04,0- (15) 

We also write 

'^j '■= (aiji "2j, a3,i, a4j) and b-, := (6ij, 62^, &3,i, ^4j) • 

Theorem 3.6. Using the notation introduced above, the second variation of J- at the VUV bilayer 
of optimal width (0) in the direction p ^Vb is given by 



de2 



£=0 



where 



Bo (flo, Sb) ■= 

iSb {-Oi - + ai_oa2,o + 03,004,0 - 4ai_o03,o + 3a2,oa3,o + 3ai,o04,o - 202,004,0} 
and, for j e N>o, 

Bj (Qji bj, duv, dvo, L) 

[duv {al^ + al^ + bi^^ + + d,o {alj + + bl^ + 6^ j }] 



— 2 {oija2,j + 03ja4j + 6i,jfo2,j + ^3,i&4,j} e 
+ 4 {oi,,03,,+5i,,63,,}e-^"''^/^ 

- 2 {01 ja4,j + 02,j03,j + &l,j64,j + &2,j&3,j} 6^^'^''''^/'^ 



+ {02,j04,j +62,jkj}e-""^''^/^ 



The proof is given in Appendix [Cl 
3.3 Variations for a monolayer 

Analogous results also hold for monolayers as defined below. In the current subsection we will 
state them. Since the proofs are completely analogous to the proofs for bilayers, we will not write 
them out here. 

The monolayer of optimal width is a structure given by functions (uq, wq) € /C with 



"0 := Xi 



.x[o,i„ 



and 



^^0 XTi,x[5„,2<5™]i 



where [T^ 



:(co + Cu + Cy). 



(16) 



(17) 



The set of admissible boundary perturbations of this structure is again restricted by regularity 
and the equal-mass constraint: 
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Definition 3.7. The set of admissible perturbations is characterised by 

V,n:={pe{W'^'{TL))' : f {p2 - Pi) ^ f (^3-^2)) 

For p G Vrn ciiT-d £ > we define a perturbed structure {u^, v^), 

, , / 1 i/x2 e (epi(a;i),(S„ +£p2(a;i)), 
[_ U otrierunse, 

, , / 1 1/.T2 € {6„, + ep2{xi),26„, + epsixi)) , 
I U otrierunse. 

We also define the subset of mass preserving perturbations: 



A picture of a monolayer of optimal width with perturbations p is shown in Figure [9] 

25m +ep3(.Ti) 

<5m + £^2(2^1) 
epi{xi) 



■m-l (P2-Pl)=/ (P3~P2)=0 



(18) 



Figure 9: The monolayer of optimal width with perturbations 

Stationarity for the monolayer of optimal width is defined analogously to stationarity for the 
bilayer, see Definition 13.31 

Lemma 3.8. The monolayer of optimal width is stationary with respect to all p E ■ 




Proof. Analogous to the proof of Lemma l3.4l we find that the first variation of the interfaces with 
respect to all p € Vm is zero. 

For any rj G C(]R) we compute 



de 



de 



rj{x2)u^{x) dx 



Sl 



= rii^rn) I P2- '7(0) / Pi, 



ri{x2)ve (x) dx 



Sl 



e=0 



■n{2Sjn) I Ps- Vi^m) / P2- 



With G the Green's function of —A on Sl from Theorem 15.11 we make the choice rj{x2) = 
G * {uq — vq){xi, X2), which is independent of xi. Using ri{Sm) ~ and 77(0) = — 77(2(5,„) > 0, we 
compute, as in the above mentioned proof, 



de 







(19) 



□ 



Note that by equation (jl9p the monolayer of optimal width is not stable with respect to 
perturbations that are allowed to change the total mass, i.e with respect to p € Vm \ Vm ■ 
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Remark 3.9. In complete analogy to Remark 13.51 we see that the optimal width condition 
(jl7p plays no role in the stationarity of the monolayer under the functional J- ^ but it plays an 
important role when considering the stationarity of the monolayer under T jM.^ the energy per 
unit mass functional. Appendix \^ also argues the relevance of optimal width when considering 
energy per unit mass. Using, |17| . 



J^(uo, Wo) = 2i:(co + c„ + c„) + ^L(5^j, 



a computation analogous to that in Remark 13.51 shows that the optimal width condition p7p is 
not necessary for stationarity under T j M with respect to the mass preserving perturbations in 
, but it is for stationarity with respect to perturbations in Vm \ T-'m- 

Similar to (|14p we express a p € Vm in terms of its Fourier modes j and 6^ j- and introduce 
the notation 

Qj := (aij, a2j,a3j) and b-, := (&i ^ , 62J, fesj) • 

Theorem 3.10. Using the notation given above, the second variation of J- at {uo,Vq) in the 
direction of p G Vm is given by 

~ = Mq {ap, 6,n) + Mj (oj, bj,duo, duv,dvo,L) , 



where 

and, for j £ N, 



(ao, 5m) ■= Sm (ai,o - a3,o)^ , 



Mj {ttj , bj , duo , duv 7 dvo , L) 



L2 

L 

+ — 



[duO {a? J + bl j} + duv {a\j + bl j] + dvo {al j + 63 



2(1 



L 



- 2 {aija2j + a2,ja3,j + bi,jb2,j + b2,jb3^j} e '^'^^"-^/^ 
Proof. Analogous to the proof of Theorem 13.61 □ 

4 Stability 

In this section we study stability of monolayers and bilayers with respect to the admissible per- 
turbations. The bilayer will be treated in Section 321 the monolayer in Section [ 



4.1 Preliminary definitions and results 

In this paper we only consider linear stability — whenever we use the words stable or unstable, this 
refers to the sign of the second derivative: 

Definition 4.1. Using the notation of Section\^ the VUV bilayer (monolayer) of optimal width 
{u(),vo) is called stable iff 

^nUe,Ve) >0, 

for every p e V^^ CPm )' ^'^'^ unstable otherwise. 
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The following property simplifies the study of stability of the bilayers and monolayers. 
Lemma 4.2. Using the notation from Theorem \ 3. 6] we have, for any y G K.** and for j > 1, 

Bj (x, y, duv,dyQ,L) = Bi (x, y, duv,dyQ,L/i) , 
Bj (x, 0, duv,dvo,L) = Bj (0, x, duv,dvo, L) , 

Bj {x,y,duv,dyo,L) = Bj (x, 0, fi„o, L) + Bj {0,y,duv,dvo, L) . 

Similarly, in the notation from Theorem \3.10\ we have, for any a;, j/ G R'^ and for j > 1, 

Mj (x, y, duo, duv,dyo,L) = Mi {x, y, duO: rf™, dyo,L/j) , 
Mj {x, 0, duQ, duv, dvo, L) ~ Mj (0, x, duo, duv, d^a, L) , 

Mj {x, y, duo, duv, d^Q, L) — Mj {x, 0, duo, duv, dva, L) + Mj (0, y, c?„o, d.^^, dyo, L) . 

Proof. These properties follow from the definitions of Bj in Theorem 13.61 and Mj in Theorem 13. 101 

□ 

4.2 Stability of the bilayer 

Throughout this subsection we will use the notation as introduced in Section 13.21 Lemma 14.21 
provides us with a simpler characterisation of stability: 

Corollary 4.3. The VUV bilayer is stable iff 

1. -Bo(flOj'^fc) ^ for all Oq G M."^ satisfying \15\). and 

2. Bi{x, 0, duv, dvo, L/j) > for all x and all j > 1. 

We therefore study Bq and Bi as quadratic forms on R"' subject to ([T5|) and investigate their 
sign. Note that Bq and Bi can be identified with symmetric 4x4 matrices, and we will continuously 
make this identification. Among other things that means we can speak of eigenvalues of Bq and 
Bi, and relate the sign of the quadratic forms to the signs of their eigenvalues. 

Lemma 4.4. Bo{a, Si,) > for all Sb > and for all Oo G R'' satisfying HIS]) . 

Proof. The Lemma follows immediately from writing Bq as 

1 1 2 1 2 1 2 

— Bo(ao,i5&) = -7: (2ai,o - 02,0 + 2a3,o - 04,0) +77(01,0-^2,0-03,0 + 04,0) +0(01,0 + 03,0) ■ 
4dfc 2 2 2 

□ 

Lemma 4.5. Two of the four eigenvalues of Bi are nonnegative for all duv, dvQ, and L; the other 
two do not have a definite sign. Denote the smallest eigenvalue by X\{duv, dvo, L) . Define 



^ _ duv _ Cu + Cv , . 



duv + dvQ Co + c„ + 2c„ 
There exists a function Ci G C([0, 1]) (see i23\) ) such that 

X\iduv,dvo,L) > C > Ci(w)- 



Proof. Note that w G (0, 1) and, by conditions (O, C G 
Let x G R"'. We now write 



2 2(co+c„+2c„) ' 2 2(co+Cu+2c„ 



C [0,1] 



2L 

Bi (x, 0, duv,dvo, L) = — Bi (x, C, v) 

TT 
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where 



^1 (x, c, log' V (c {xl + xl) + (1 - c) {4 + 4)) (21) 



+ (l + logw) H + a:^) + 1 + 2,2^ 



(a;ia;2 + X3X4) v + 2xiXzv^ — (2:1X4 + X2X3) + -X2Xiv'^. 



Note that when xi ~ X'i ~ 0, 



1 '2 , 1 



Bi {X, C, = (1 - C) (4 + ^4) + 4 (a^2 + xl) + -X2XiV^ > 0, 

so that by the max-min characterisation of the third eigenvalue A3, for fixed (^,v, we have 
A3 = max min Bi{x,(^,v) > min C, > 0, 

dimL=2 xeW^/L xi=X3=0 
\x\ = l 1^1 = 1 

implying that the largest two eigenvalues are always non-negative. 

We now turn to the question of existence of admissible x such that Bi is negative, and we 
simplify the problem by minimizing Bi with respect to X2 and X4 under fixed xi and 0:3. The 
stationarity conditions gf^-Si (x, (,v) = and gf^-Bi (a;, C,v) = lead to the equations 

xT ) ^ detA(C,«)^^^'"^ « j U3 

where 

i-^(l-C)log3i^ -^v^ 



Inserting these results into Bi gives 



Si {xi,X2^\x3,xl^\(,v) = {xi,X3) B{C,v) {X1,X3) , 



where the matrix entries of B are given by 

^ * 1 

= 522(C,w) = log -Clog u 



Sl2(C,w) =B2l(C,«) 



5c 

The eigenvalues of B are 



(3(-l + v^) - 4(-l + C) log^ vj (3(-l + u^) - 4(-l + C) log^ 
9(-l + + 8(-l + C) (-3 - 2(-l + C) log^ w) log^ V 
(3v{-l + v^)- 4tj(-l + C) log^ ^ 
9(-l + + 8(-l + C) (-3 - 2(-l + C) log' log' 



G_(C,z;) :=l-z;^+log^;~^Clog'z;+ ^ , ^^'^ 3 



3 ° 3(-l + t;4)-4(-l + C)log''u 
= (3(-l + v^) - 4(-l + C) log=^ v) "\„(C, ^;), 

2 , l.,„„3.. 3^.2(1 + ^,2)2 



G+(C,w) := l+u^ + logu- -Clog I. 



3 3(l + u4) +4(-l + C)log'^u 

= (3(1 + + 4(-l + C) V "'/^+(C, ^^), 
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with 

h.{C,v) := (^log'w) C'+ (-^ log' w-41og%+ (-3 + 4^2 -^;4)log3^;^C 

- 3(1 - v^)^ + 3(-l + V*) logv + 4(1 - v^) log^ V + 41og^ v, 
h+iC, v) := - Q log' C' + Q log' V + 41og^ u + (3 + 4^2 _ ^;4) j^gS ^ 
+ 3(1 - v^) + 3(1 + v"^) logv - 4(1 + v^) log^ u - 4 log"^ w. 
Note that G_ < G+, smce for v € (0, 1), ( G [0, 1], 

3(1 + v^) + 4(-l + C) log^ f > 0, 3(-l + v^) - 4(-l + C) log^ f < 0, 

and thus 

-2i;2(3(-l + -4(-l + log^u)^ 

G+iCv) - G_(C, v) = ^ 1 > 0. 

+ w4) + 4(-l + C) log^ vj + v^) - 4(-l + C) log^ vj 

We have now the foUowing equivalences: 

Vx G M*,Bi {x,0,duv,dyo,L) > 0^=^\fx e R'^,Bi (x,C,u) > 

{C,v)>0. 

We prove the following characterisation of the sign of G_ : 

G_(C,i;) >O^C>Ci(^), (22) 

where 

Ci(f) = (81og^u)"i^9- 12t;2 + 3u'^+ (41ogt;)(3 + log2u) 

+ 1 225 - 504^2 + 342^/ - 72v^ + 9v^ + (360 - 288^^ - 72v*) logv 

+144 log^ V + (-120 + 96u^ + 24t/) log^ u - 96 log'^ w + 16 log^ ^} ' ^ ■ 

(23) 

The details of this calculation can be found in Appendix |D] This concludes the proof. □ 
The function gi mentioned in Theorem 11.21 in the introduction is related to , given in (|23p , 



by 

9^{e):=C^{e'^/'). (24) 



Remark 4.6. The four eigenvalues of Bi from the proof of Lemma are 
^ ^45 - 36u2 - 9^2 + 36 log u - 12 log^ v 

± I (^-45 + 36v^ + 9v^ - 36 log u + 12 log^ ^ 

- 144(9- 18i;2 + 9w"* + 91og-!;-9u"*logi;- 12 log^ u + 12-!;2 log^ w + 9Clog^ u 

- 12v^C log^ V + 3v*C log^ w - 12 log u"* + 12c log'^ v + 4C log^ v - 4^ log' v 
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0.6 



1^ 



G+ > 




0.2 0.4 0.6 



(a) The sign of G_ 



0.2 0.4 0.6 0.8 

(b) The sign of G+ 



Figurc 10: The sign in parameter space of the eigenvalues G_ < G+. The boundary between the 
two regions in the left-hand figure is given by C = ^ (u) . 



and 
72 I 



45 + mv^ + 9v^ + 36 log w - 12 log'^ v 



±{(45 + 36t;2 + 9u'' + 36 log u - 12 log^ vj 

144 (^-9 + 9u** - 9 log u - 9u^ log u + 12 log^ v + 12v^ log^ v - 9C log^ v 

- l2v^C log^ V + 3v^C log^ f + 12 log - 12C log'' v - AC log^ v + 4C^ log^ v) 



Plotting the areas where these eigenvalues are negative shows that the eigenvalues with the plus 
sign chosen for ± are positive everywhere for u S (0, 1) and C G [Oi !]• The plots for the other two 
eigenvalues correspond with those in Figure 1101 

Collecting Lemmas 14.41 and 14.51 we can summarise the stability properties with the use of 
Corollarv l4.3l as follows: 



Theorem 4.7. Let v, and Ci be as in Lemma \4-5\ Define the functions C^.j J > 1, and C. by 

Then the VUV bilayer of optimal width 0j is stable with respect to all (mass-conserving) pertur- 
bations in V^^ iff 

C > Civ)- 



This is Theorem 11.21 from the introduction. Its implications are illustrated in Figure I4b[ 

Remark 4.8. Note that the statement in Thcorcm l4. 71 about the positivity of the second variation 
also holds true if we allow the perturbations to come from the larger set of perturbations Vb, instead 
of V^^ ■ However, as stated in Remark 13.21 the bilayer of optimal width is not stationary under 
perturbations that do not preserve mass. 

We next show that ( is bounded from above away from 1. Therefore there is a threshold a (as 
mentioned in the introduction) such that the bilayer is stable if C > a. 
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Lemma 4.9. Let C he as in Theorem \4.7\ then there exists a G (0, 1) such that for all v G (0, 1), 

C(w) < a < 1. 

Proof. First note that per definition of ( it suffices to show that there exists a c G (0, 1), such that 
for all uG (0,1), 

Ci(w) < c < 1. 

Since (i is continuous on the interval (0, 1) and goes to zero for J, and to | — \ for u t 1 
(see Remark Id. ip . this is equivalent to 

(81og3i;)(Ci(i;)-l) >0. 

By (gOl) we know that 



0< I (^9-12?j2 + 3w'' + (41ogu)(3 + log^j)) -81of 



< 225 - 504^2 + M2v^ - 72^6 + 9i;** + (360 - 288^2 - 72v^) logv 

+ 144 log^ V + (-120 + 96-u2 + 241;"*) log^ v - 961og^ v + 161og^ v. 

Taking square roots completes the proof. □ 

Remark 4.10. To find out the stable and unstable first-order Fourier modes of deformation 
for the bilayer, we compute the eigenvectors belonging to the positive and (potentially) negative 
eigenvalues of Bi from ([21]) . For the stable directions we find 



\12u(l + u2)V J V f^{C^v) + ^fliCv) I 



ar (C, V) ^^-4-^ (.i(C, V) - v/^) , -1, "i + ff^,! , 

yi2v{-l + v^)\ I V gsiCv)- ^g4C^ J 

where 

/i (C, u) := -9 - 12^2 + 31;"* - 12 log u + (-4 + 8C) log^ v 
/2(C, v) := -9 - 3^2(6 + v^)- 12 log u + (-4 + 8C) log^ v 
/3(C,w) := 15 + 3^2(4 + ,7^) - 121ogw + (-4 + 8C)log3w 
/4(C, v) := 9 (9 + 40^2 + 42^4 + 8u^ + v^) 

gv(-3 + (-1 + 20 log^ v) (3(-3 - 4u^ + v^) - 61ogt; + (-2 + 4C) log^ v 



gi{C,v) := -9 + Uv^ - Sv"^ - 12 log w + (-4 + 8C) log^ v 
32 (C, v):=9- 3v^{2 + w^) + 12 logu + (4 - 8C) log^ v 
gaiC, v) := -15 + 3^2(4 + v^) + 12 logv + (4 - 8C) log^ v 
giiC, v) := 9(-l + i>2)2(i + ^,2)(9 ^ ^2) 

+ 8 logu(-3 + (-1 + 2C) log^ (-3(3 - Av^ + v^) - 6 logv + (-2 + 4C) log^ . 

The directions belonging to the eigenvalues that can become negative, corresponding to the eigen- 

* i-r-FH 

values G+ and G- of the reduced matrix B in the proof of Lemma 14.51 are 



^) - ( iMTT^ i^^^^' + ' ^ /3(C^)-V/.(C,.) ' V ' 

\ 12i;(-l + i;2) V y "53(C,w) + \/ff4(C^ , 
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respectively. 

The direction of the perturbation a"^ is depicted in Figure [7al Here we have chosen the 
values duv = 0.7, d„o = 0.3, L = 5 and e — 0.25. Similarly we get Figures Ffbl [7cl and [7dl using 
perturbations a"^, and a^^ respectively. 



4.3 Stability of the monolayer 

We now redo the arguments for the monolayer of optimal width (jl6p . Throughout this subsection 
we use the notation of Section 13.31 

We can simplify Mi a bit by writing 

^ g-27r(5„/L '^"0 _ C„ + Co dyo _ Ci, + Cq 

duO + duv + dyo 2(co + + dyo 2(co + c„ + 

Note the slightly different definition of i/ than for the bilayer ([20|) . Then, for all a; e R'^, 

Ml {x,0,duo,duv,dvo,L) = —Mi {x, q,<,, u) , 

TT 



where 



We now can write 
with 
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Ml [x, Q, ^, v) := log^ {qxI + Q- <C)xl + <ra;^) 
+ 2(l + logi.)a;^ + i(a;2+x2) 
— 2 (xi + + 0:2X3) V + xixs^/^. 

Ml (x, e, ^, I/) = x'^M{g, ^, I/) X, 



-|elog^:^ + i ~iy 



2 2 

M(£), ^, I/) ( -iy -|(l-g-cr)log'''i/ + 2(l + logz/) -i^ 

ii^^ -i^ -|^log^j^ + i 

This matrix is well defined for all g,(; E M., > 0, but note that the nonnegativity of the parameters 
Co, Ck, and Cy — or equivalently conditions (O — translates into 

and furthermore v € (0, 1) by definition. 

Remark 4.11. As explained in the introduction and Appendix [B] we assume throughout 
the paper that for the monolayer the interfaces U-0 and V-0 are penalised equally strongly, i.e. 
duO = dyo or equivalently c„ = Cy. Under this assumption x = k, and the inequalities in (|25p 
imply that x ^-nd k take values in [|:, 

Lemma 4.12. Let Cy = Cy. Two of the three eigenvalues of M v) are nonnegative for all 
V E (0, 1) and E [i^h]- The third eigenvalue is given by 



^ {ei{<;,v) - ^62(^,1/)) , 



where v £ (0, 1), ^ £ [j, ^] and ei and 62 are given in ^27\ ) and i28\). The sign of is characterised 
by 

E^{<,,y)>Q^^<^2{v) (26) 

with (,2 as given in II 29]). 
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Proof. Since we are interested in the case where Cu = Cv we will take p = "7 from here on. which 
turns the conditions into | < ^ < ^. For the three eigenvalues of Mi((^,i;,u) we compute 



6 



1.3 (t, i^) ^ (ei(<;, ly) ± ^/e^{syj^ , 
where 

ei(<?,t^) := 15 + 3i^2 _^ (12 - 4 log^ + 4^ log^ t^) log z^, (27) 
e2('? , i^) := 81 + 2341^2 + 9i,4 ^ 2I6 log - 72i^2 ^ _,_ ^44 log^ 
- 72 log^ 1/ + 24i/^ log^ - 96 log"* + 16 log^ 1^ 
+ (216 log^ V - 72i/2 log^ zy + 288 log"^ ly - 96 log^ v) ^ 
+ (l441og<'z.)<;2 (28) 

and we choose the plus sign for E2 and the minus sign for E3. 

First note that i' £ (0, 1) and x ^ imply that Ei is always positive. i?2,3 are real, since they 
are the eigenvalues of a symmetric matrix and thus e2(<j, v) > for all S M and for all ly S (0, 1). 

Since for all x > and x < 1/2 we have (1 - x)x^ ~ 3x > {l/2)x^ - 3x > -2^2, 

ei{<;, = + + 4[(1 - x)l logz^l^ - 3| logz^l] > 15 - 8^/2 > 0. 

Combining this result with e2(<j,i^) > 0, we conclude that £'2(1;, i^) > for all admissible <r, i/. 
Thus, the only eigenvalue that might be negative (in all or part of parameter space) is E3. 
To prove the statements in (j26p we compute 

e2(<r, ly)) = 9(1 - ly^) + 9(1 + ly'^) log - 3(1 + ly'^) log^ ly 
+ (6(-l + ly^) log3 ly + 4(-3 + log2 ly) log'' ly^ 
- (Slog^i.) ^2. 

This expression is negative on (0, 1) if and only if <; G [j, <ri(i')) U ('?2('^)7 5] and zero if and only 
if IT = <;i{iy) or <j = <?2('^), where 

TTT^ (^('^) ± , (29) 

16 log ly ^ ' 



with 



!{v) (6(-l -^v-)+ 4(-3 + log^ v) log i^) log^ i^; 

g{v) 96 log*^ 1/(^3(1 - j/^) + 3(1 + i/^) log - (1 + i^^) log^ 

+ (6(-l + v^) Xoi V + 4(log^ v-i) log" 

The minus sign is chosen in while in c:2 we choose the plus sign. Plots of "^i and ^2 are shown in 
Figure [11] 

It is left to prove now that ^1(1^) < 1/4 for all v G (0, 1). We will actually prove the stronger 
statement q(i^) < 0, which follows from 

g[v) > for < J/ < 1, 

^ !{vf - 3{v) < for < 1/ < 1, 

<=4> 3(l-i^2) + 3(l + i/2)logj/-(l + t^^)log^t^>0 forO<j/<l, 

l-z/2 

<^=^ 3- J + 3 log J/ - log-^ i^>0 forO<j/<l, 



W — — log ly Q 

<^=J> 3 tanh w — 3u' + ui > for w > 0. 
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To prove that this last inequahty holds, wc define h{w) := tanhu; — 3w + and use tanh' w = 
1 — tanh^w, to compute that h"'{w) = Gtanh^w (—3 tanh'' w + 4) > 0. From this it follows by 
integration that h{w) > for all w > 0. □ 



?2 



0.4 0.6 O.E 



0.4 
0.3 
0.2 



(a) Plot of fi away from v = 1 



(b) Plot of f2, with the admissible range 

?2 e [ ■ 



indicated 



Figure 11: 



Remark 4.13. For the excluded endpoints and 1 we find 



limcTi = 0, 
lim^2 = ^, 



lim<Ji = — o 
hmi;2 = -■ 



The limits for v '\ 1 were found by calculating the first terms in the Taylor expansion of q^2- 

Figure fT2l shows the parts of parameter space where is positive and negative, both on the 
'I 1^ 



admissible domain f^, for <j as well as extended to (0, 1) 




0.2 0.4 0.6 O.E 
(a) Sign of E3 



E3<0 



U U.^ U.4 U.b O.d 1 

(b) Sign of i?3 extended to (0, 1) for ? 



Figure 12: 



Remark 4.14. Expanding E3 around v = 1 gives 

4 



E,{,, u) = -(1 - 50(1 -yf+O ((1 - yf) 



ioi V ] \. Since 1- 5^ < -i for e we can conclude that for v close to 1 (or equivalently 

large L) the monolayer is unstable for all admissible values for the interfacial (surface tensions) 



coefficients dij (or q). This corresponds to what is shown in Figure n-2al 

Taking into account the assumption d^Q = d„o, the condition 1 — 5<; < for negativity of E^ is 
equivalent to duv < ^{duo + dvo). In [171 Theorem 8] we show that for a circular two-dimensional 
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monolayer the term in in the energy per unit mass T jM that is quadratic in the curvature is given 
by 

1 



where m is the thickness of the layers and k is the curvature. Taking m = (5m we find that this 
term becomes negative exactly as duv < |(c^uO + d^o), showing that the (large) circular monolayer 
loses stability at the same point as the flat monolayer on large domains. Note that conditions ^ 
imply duv < |(c^«o + dvo)- 

In order to compare the results for the monolayer to those for the bilayer, we introduce the 
relative U-V interface penalisation 



1 - £» - 



2(co + Cu + Cv) ' 



(30) 



analogous to C for the bilayer in Lemma [1751 Note that conditions ^ give yU £ [0, i]. In terms of 
the surface tension coefficients, 

duv 

duO + duv + dvo ' 
fi is interpreted as the relative penalisation of the U-V interface. 

Theorem 4.15. Let duo — dvo and let ^2 be as in Lemma \4.1^ Define the functions and <f by 

S.Ai^) ■= '^2(1^^), 'f:=inf£., fi:=l-2q. 
J j>i ■' 

The monolayer of optimal width il6]) is stable with respect to perturbations in if and only if 

Proof. First we work with <j as in Lemma 14.121 and afterwards we translate the results into con- 
ditions on ju. By Definition 14.11 and Theorem 13.101 in order to prove stability, we have to prove 
that 

00 

Mo (ao, 5m) + ^ Mj (ttj, bj, d„o, duv,dvo,L) > 0, 

i=i 

for all admissible perturbations. Per definition we have Mq (oo,^m) := Sm (ai.o — tJs.o)^ > 0. By 
Lemma [4.121 we know that if (1^,1;) is such that c; e [^,^2(1^)] then AIi {ai,0,duo, duv, dvo,L) > 
for all p G V,n- By Lemma we have for j > I, 

Vaj, bj, Mj (oj, bj,duo, duv,dvo, L) > <;=^ Vai, Mi (ai, 0, d L/j) > 0, 

thus we see that, if <; e [■|,if(;^)] is satisfied, then for all j > 1, for all a^, and for all bj, 
Mj (qj, bj, duo, duv, dvo, L) > 0. 

Now note that ? = and thus 



which proves the statement of the theorem. □ 
To make the connection to Theorem 11.11 in the introduction, the function /i is defined by 



(31) 



where X2 is given in (P^ . 

Figure [4a| illustrates the stability properties of the monolayer of optimal width from Theo- 
rem UTTHl 
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Remark 4.16. In Theorem 14.151 wc only consider perturbations in 7^,^/, i.e. perturbations that 
keep the total mass fixed. The statement about the positivity of the second variation stiU holds 
if we consider the larger set of perturbations Vm, however, for these perturbations the monolayer 
of optimal width is not a stationary point, as was noted after Lemma 13.81 

Remark 4.17. To find the stable and unstable first order Fourier modes of deformation we com- 
pute the eigenvectors belonging to the positive eigenvalues of Mi (ai, i;, ^, i^) and to the eigenvalues 
that are negative for some parameter choices. For the positive, stable directions we find 

a\'{^,y) := (-1,0,1), 

ar(^,i.) := (^1, J-(/ii(^,z.) + y)I^(^),l^, 

where 

hi{q,v) := -'i + ?,v'^ -\2 log v+{A-l2q) log^ v, 
h2{<;,v) 9 (9 + 26i/2 + i/'i) 

+ 8 log (3 + (-1 + 3<r) log^ v) (9 - 3i^2 + 6 log v+{-2 + 6(^) log^ v) . 

The direction belonging to the eigenvalues that can become negative, corresponding to the eigen- 
value i?3 of in Lemma I4.12( is 

Figure I5bl shows the monolayer with a perturbation corresponding to a^^ . Here we have chosen 
the values duo = 1, duv — 0.7, d^o ~ 0.3, L ~ 5 and e = 0.25. Similarly we get Figure [5c] using 
a^^, and Figure [5al using a". 



4.4 Discussion and comparison 



In Sections l4.2l and l4.3l we found conditions for the stability of monolayers and bilayers with respect 
to some admissible perturbations. The main results are visualised in Figures [13] and [14] for the 
monolayer and Figure [TSl for the bilayer. 




5 10 15 20 25 30 35 

(a) Plot of fi as a function of L j Sm 
27r(log!/-l)-l 



L/S„ 




5 10 15 20 25 30 35 
(b) Plot of fi restricted to the admissible 
range fi e [0, |] 



Figure 13: For the plots of /i from Theorem 14. 151 we have approximated <f by mini<j<ioo£- 



The monolayer is stable with respect to perturbations of the interface ii fJ.> p. f Theorem 14. 15p . 
and the bilayer is stable with respect to mass-preserving perturbations of the interface if C ^ C 
(Theorem ]12]). 

fl and ^ display very similar overall behaviour. They both rapidly increase for small values of 
L/Sm or L/Sb until they settle down around a value for /it or ^ close to 0.6. Around this value both 
p, and C oscillate as with increasing L different Fourier modes become dominant. The similarity 
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0.6098 
0.6096 
0.6094 
0.6092 
0.609 
0.6088 



Figure 14: Plot of /2 from Theorem 14.151 as a function of L/5m = 27r(logi^ ^) ^, showing the 
smaU-scalc oscillations where different Fourier orders become the dominant contributors. 




c 



(a) Plot of as a function of L/Si, 
27r{logu~i)~l 




L/Sb 



0.6569 
0.6558 
0.6557 
0.6556 
0.6555 

(b) Plot of f showing the small scale oscil- 
lations where different Fourier orders be- 
come the dominant contributors 



Figure 15: For the plots of ( from Theorem 14.71 we have approximated C by maxi<j<ioo C 



is broken, however, by the restriction of /i to [O, ^] . Because of this the monolayer is unstable for 
all values of L/6m greater than about 6 (see Figure [T3b)) . while the bilayer can be stable for all 
values of L/St (see Figure [T5a|) . 

Remark that higher relative penalisation of the U-Vinterfaces, i.e. higher values of /i and 
improves stability. For the bilayer a sufhciently high value of C even guarantees stability in the 
sense discussed here (Lemma 14. 9p . This reinforces the notion that duv pla-ys a special role in the 
diblock copolymer- homopolymer problem, which was also encountered also in [17|, §2.2]. 

5 Green's function on a periodic two-dimensional strip 

When computing the first and second variation of for monolayers and bilayers in Section [3] 
we required an explicit formula for the Green's function of —A on Sl- We present this Green's 
function here and prove that it satisfies the necessary conditions. For a heuristic derivation we 
refer to [Ml §6.3.1]. 

Theorem 5.1. Define G : Sl\ {(0,0)} R in L]^^{Sl) as follows: 

G{xi,X2) := ^ log (^2 cosh (^^^ _2cos(^^^^ . (32) 

Then the equation ~AG{xi, X2) = d{xi, X2) is satisfied with periodic boundary conditions G{0, X2) = 
G{L,X2) and ■^^G{0, X2) = -^G{L, X2). Writing the Fourier expansion of G in xi gives 

q—1 ^ 
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Proof. Wc first prove that G, as given in equation ([5^ . satisfies the equation — AG(.ti, 2:2) ~ S{xi, X2) 
in tlie sense of distributions, i.e. we show that for all cj) £ C^{Sl), 



G{xi,X2)i-A(t>){xi,X2) dxidx2 = 0(0,0). 



Sl 



Note that the constant term — ^ log 2 implicitly present in ([5^ as the factor 2 in the logarithm 
is of no importance here and so we will leave it out of subsequent calculations 0. 
We write 



J Sr. 



'Sl 

lim 



lim 



G{-A<p) dC^ 



Sl\B{0.e} 



i- f GV<t)-vd'H^- f 

\ JdB{0,e) Js 



AG4)dC^ + 



SL\B(0,e) 



VG-v4)dn^ , 



dB(0,e 



where 5(0, e) is the closed ball of radius e and with the origin as center, v is the unit outward 
normal to Sl \ -6(0, e), which means v points into _B(0,e). Denote the three terms by Je, and 
A'e respectively. The integral vanishes: 



lim 141 < lim ||V(/)||l- / \G\ dH^ 

JdB{0,e} 



lim ||V0||loc, 27r£ 



log(^(e2 + 0(e^)) 



= 0. 



For Jg we calculate 
yG{xi,x2) 



1 

2L 



cosh 



2t:x2 



ilTXl 

L 



sin(2£ii) 
sinh(^) 



For notational convenience we will write C{xx,X2) '■— cosh(2^S3.) — cos ( ^^^ )■ Then we can 
compute that at (a;i,X2) ^ (0,0) 



X 

-^G{xi,X2) 

oxi, 



^^^G{X^,X2)- ^2 



l2 



TT 
I2 



C(a;i,a;2) ^ sin^ 
C(a;i, X2)~'^ sinh^ 



27ra;i 
/ 27ra;2 



C(a;i,a;2) ^ cos 



/ 27rxi 



C{x\,X2) ^ cosh 



/ 2tTX2 



which gives AG{xi, X2) = 0, from which it follows that 4 = for all e > 0. 

To determine lime_,o we approximate G by Gr2 (xi, X2) = — (47r)^^ log(xi +2:2)1 the Green's 
function of —A on R^. Estimating the difference on dB{Q,e) by 



|VG(xi,X2) - VGR2(a;i, 0:2)1 



we calculate 



1 ^lI^g^ + ^g^ + o{{xl+xir/^) ^ 1 xiei+X2e2 



2L 2^(xl+xl)+0{{xl + xiy) 2tt xj + x 
0{{xl + xl)^/'^) as xl + xl = s'^ 0, 



lim = lim / V(G - Gr2 ) • i^cf) dV} + hm / VGr2 • dH^ 

e^O ^^OJdB{a,e) J dB{0,e) 

1 



= + hm 



xiei+X2e2 xiei + a;2e2 



(xi, 2:2) d7i^(a;i, 2:2) = 0(0, 0). 



*The reason for adding it in I I32I I in the first place is to get a Fourier series without a term independent of xi 
and X2- 
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Taking these results together shows that \im^io{I^ + + K^) — 0(0, 0) and thus —AG = 5 holds 
in the sense of distributions. 

To prove that the Fourier series in ([55]) corresponds to the Green's function (|32p . let G be 
given by ([5^ and G by ([551) . Note that for every X2 7^ the series converges absolutely: 



E 

9=1 



1 

-e 
9 



If a;2 7^ we then calculate 

-e ^ cos 
9=1 ^ 



2'KXiq 



^ 1 



e 



9=1 



-|2;2|+ia;i 



1-e 



-Relog(l-eT(-l^2|+«xi)^ 
-logll - eT(-l^2|+ia;i)| 



- - log I 1 — 2e i cos 



2'KXl 



= --log(2e- 



1 2^\^2\ 



1 



2^2 



log 2 cosh 



27r|x2| 
L 



- e i 
27ra;i 

2 cos 



1 



27ra;i 



L 



(34) 



so that the partial sums X]g=i cos i^ ^^"^^ ) converge pointwisc to 27rG'(xi, X2) + ■j^\x2\ 

for almost all (xi,X2) ^ Sl- Since the partial sums are all bounded by the Lj^Qj.-function on the 
right hand side of ([M)) the Dominated Convergence Theorem yields G S L'^{Sl)- Together with 
G = G a.c. on Sl this shows that G = G in ^^^^(5^). □ 

Corollary 5.2. Let G be as in ^ and let X2eR\{0}. Then 



Proof. For all q > 1, 



G(a;i, 0:2) dxi 
f 27rga;i 



cos 



V L 



dxi = 0. 



□ 



Also note that G(— xi, 2:2) = G(xi, X2) and G(xi, —X2) = G{xi, X2)- 



6 Discussion and conclusions 
6.1 Comparing mono- and bilayers 

In this paper we showed that bilayers can be both stable and unstable, depending on the pa- 
rameters: when the U-V interface penalty is strong enough relative to the penalties of the other 
interfaces, the bilayer is stable. On the other hand, monolayers are unstable as soon as the strip 
is wide enough to accommodate the unstable wavelengths, regardless of the values of the interface 
penalisation. 

The bilayer can be thought of as two juxtaposed monolayers, and therefore the question presents 
itself how the unstable mode of the monolayer is prevented in the bilayer context. The correct 
answer seems to be that the unstable mode is actually not prevented at all; it continues to exist 
in the context of the bilayer, as can be witnessed in Figures [5al and (especially) [Tbl 
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The reason why this unstable mode docs not make every bilayer unstable lies in the admissible 
values of the coefBcicnts, which are different in the two cases. For the VUV bilayer, for instance, 
the value of the U-0 interface penalty duo is irrelevant; therefore, by choosing d^o '■= dyo + duv, 
every choice of d^v and d„o becomes admissible, and most importantly, the case of purely U-V 
penalisation (C w 1, or d^o « 0) is therefore allowed. For the monolayer, however, the conditions ([3]) 
imply that the two side interfaces (0-U and V-0) are necessarily penalised at least half as strongly as 
the central (U-V) interface. Most of the white (stable) region in Figure [T2bl therefore is inaccessible, 
and only the unstable region remains as can be seen in Figure fl2al (Figures 1121 onlv show stability 
of the first Fourier mode, but the situation is similar for the higher modes). 



6.2 Comparison with |34j 

In previous work [34] one of the authors (Pelctier) and Roger studied a related functional, 

e I |Vu| H — di{u,v) if(M, w)£/Ce, 
ge{u,v) := { Jr^ ^ (35) 

xj otherwise. 
Here di{-, •) is the Monge-Kantorovich distance with cost function c{x,y) = \x — y\, [56| . and 

/Ce := |(u,w) e BV(M^;{0,l/e}) X Li(R2;{0,l/e}) : WW = a.c., and [ u^f v = m\ . 

Apart from the choices Cg = Cy = and c„ = 1, the main difference between and (|35p is the 
different non-local term. 

The scaling (constant mass but increasing amplitude 1/e) implies that the supports of u and v 
shrink to zero measure. The main goal in [34] was to investigate the limit e and characterise 
the limiting structures and their energy. 

The main result, a F-convergence theorem, can be interpreted as stating — in a very weak 
sense — that the limiting structures are VUV bilayers; in the limit e ^ these bilayers have a 
thickness equal to 4e and their curvature is bounded in L^. Most importantly, in connection 
with the present paper, the limit energy depends on the curvature in a stable way: the energy is 
minimal for straight bilayers and increases with curvature. 

This result compares well with the results of this paper. The functional of [M] penalises 
only U-V and U-0 interfaces; the V-0 interface is free, or in terms of this paper C = 1- Both in [M] 
and in the present paper we therefore find that bilayers of optimal width are stable, although the 
precise results and their methods of proof are very different. 



6.3 Comparison with 'wriggled lamellar' solutions 

In a series of papers [30l[37l|40] Muratov and Ren & Wei investigate the stability of one-dimensional 
layered (lamellar) structures for copolymer melts — the case u + v = 1. They find that for a critical 
value of the lamellar spacing the straight lamellar structures become unstable and a stable branch 
of curved, 'wriggled' lamellar structures bifurcates. Muratov considers unbounded domains and 
finds that the loss of stability happens at exactly the optimal value of the width: for any larger value 
of the width unstable directions exist with very large wavelength. Ren and Wei consider bounded 
domains, which provides a natural limit on the wavelength of perturbations, and consequently 
they find that at the optimal width the straight lamellar structures are stable, and the bifurcation 
occurs at slightly larger width. 

The system studied in this paper is different in that there are three types of interfaces, not 
one; for comparison purposes one can identify the pure-melt case described above with the case 
of pure U-V interface penalisation for bilayers (C = 1). In this case the bilayer of optimal width 
is stable, and this result mirrors the stability result of Ren and Wei for optimal-width lamellar 
structures. 
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6.4 Generalizations and extensions 



One might wonder whether the functional T depends in a smooth manner on the perturbations. 
The calculation of the second derivative of the functional in the melt case done by Choksi and 
Sternberg [10] suggests that the second derivative of T depends continuously on ly^'^-regular 
perturbations of the interfaces. In that case the functional T is of class C^, and the linear stability 
analysis of the current paper automatically implies the equivalent nonlinear stability properties. 

One can also wonder whether the class of perturbations that are considered — those described 
by functions of the variable x\ g — is not too restrictive. The class of all perturbations that 
are small in L^, for instance, also includes many perturbations with small inclusions of one phase 
in another, which are not covered here. We believe that these will generally be less advantageous, 
since the results of this chapter show that perturbations with fast oscillations are energetically 
expensive (because the layers are stable with respect to the admissible perturbations for most 
values of the surface tension coefficients if L is small). The same conclusion can be reached by a 
slightly different, heuristic argument as follows. Within the class of uniformly bounded functions 
the i7~^-norm is continuous with respect to the L^-topology, as can be seen from 

where (p solves —Aip — f. Therefore within that class of functions the iJ^^-norm is also continuous 
with respect to the area of the inclusion; for small inclusions, with a large circumference-to-area 
ratio, a possible decrease in the i/~^-norm is thus dwarfed by the increase in interfacial length 
associated with such an inclusion. 

Note that the problem has not completely been non-dimensionalised; it is possible to rescale the 
problem by the length scale L, resulting in a three-parameter problem (in the rescaled parameters 
Co, c„, and Cy). Instead we keep the length scale explicitly in the problem to illustrate the length- 
scale dependence of the stability properties. 



6.5 Diffuse interface model 



The functional is the sharp interface limit (via F-convergence) of a well-known diffuse-interface 
functional l28l[6l 



||V7.|2 + ||V«|2 + ||V(ii + v)\'' + e~^W{u, v) 



dx + — \\ u — v\ 



Here W is a triple-well potential with wells at (0,0), (1,0), and (0, 1). The coefficients duv: duo, 
and dyQ in the sharp interface limit depend on the specific form of W via 



ikl 



2inf 



VWhmi'it)\ dt:-i& C\[0, 1]; (R+)"), 7(0) - au,l{l) = m 



where = (1,0), = (0,1), and ao = (0,0). 

By the properties of F-convergence minimisers of converge to minimisers of T [12[ Corol- 
lary 7.17]. Therefore our results indicate that in the regions of their respective instability mono- 
layers and bilayers are not minimisers for JF^ for small e. 



A Relevance of energy per unit mass for partial localisation 

Throughout this paper we concentrate on layered structures with a specific width: the width that 
minimises the ratio of (one-dimensional) energy to (one-dimensional) mass. The origin for this 
choice lies in our interest in partially localised structures, as we now explain. 

Since we are interested in long thin structures, we might first ask ourself the question what 
minimisers of T on the full domain look like if we restrict the admissible functions to be 
rectangles with a fixed mass, oriented such that the long axis is parallel to the a;i-axis. 
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If the rectangle has a large aspect ratio, the structure is roughly constant in the xi-direction. 
We can interpret the rectangle then as a one-dimensional structure in the a:;2-direction, extended 
trivially in the xi-direction and cut off at a certain length, a. In |17| it is proven that for such a 
trivially extended one-dimensional structure the energy J- per mass M is approximately equal to 
the one-dimensional energy Fid per mass of the cross-section Mid: 

Put differently: although the energy depends on the structure in a nonlocal manner, for large 
mass (i.e. long rectangles) the energy is essentially equal to the one-dimensional energy of the 
cross-section times the length of the rectangle. Effects near the cut off points are less important. 

This implies that the miminiser of T in the class of rectangles with large constrained mass 
should have a thickness Mid such that J^id/Mid is minimal. Also when studying the stability of 
layered structures, it thus makes sense to concentrate on structures of optimal width, in the sense 
as described above. 

In a monolayer of optimal width the U- and V-layers both have width |17| 



3^1/3 



(co + Cu + c„) 



1/3 



2, 

while for the bilayer the thickness of the inner layer is 

2^6 := 6i/3(co + Cu + 2c„)i/3 (VUV) or 25b Q^'^{co + 2c„ + c„)^/^ (UVU). 



B Relevance of the choice Cu = c„ for monolayers 

The choice Cu = Cv for monolayers is similarly inspired by our interest in partial localisation and 
more or less forced upon us by the periodicity in the xi-dircction. If the U-0 and V-0 interfaces 
are penalised unequally, then a monolayer structure in likely will tend to curve, in order to 
reduce the length of the 'expensive' interface at the expense of the 'cheap' interface. 

When c„ 7^ Cu, therefore, a straight monolayer is not even stationary under perturbations that 
allow for curving of the whole monolayer. The setup in the context of the strip Sl disallows such 
curving over the whole length of the layer because of the periodicity in the a;i-direction. Therefore 
this instationarity is rendered invisible on Sl- However, with our interest in partial localisation 
in mind we make the choice c„ = Cy in the case of monolayers throughout this paper. 



C Proof of Theorem D 



For the interfacial terms we directly compute from pOl 



d£2 



JSl JSl JSl / 



e=0 



duv 



/ 2 , ,2 
Pi +P3 



/2 ,2 

P2 +Pi 



dx. (36) 



In order to compute jr^\\ 



:=0 



we split up the norm as follows: 



\Ue-Ve\\H-i(SL 



L rL 



"'0 



fe{xi,^l)dCl dxi, 
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where 



G{x - £,) d^2 dx2 



2(5b— ep4(a:i) — 2(5b — 6P4(5l) 



5t+epi{xi) rSb+epi{^i) 



G{x - d^2 dx2 



Sb+epi(xi) JSb+spi(.ii) 



G{x — ^) d£,2 dx2 



- 2 



Sb+epi (xi) 



-Sb-ep3(ii) 



'Si,-!^P3{xi) J -25b-epi{^i) 
2Si,+ep2{xi) r — Sb—epaiii) 

Sb+spi{xi) •! -25b-epi{(,i) 



G{x - d6 dx2 - 2 
G{x — ^) d^2 dx2 



Sb-ep3{xi) J -Sb-spsiii) 



25b+ep2{xi) rSb+epii^i) 



Sb+epi {xi) 



G{x - d^2 dx2 



(38) 



We compute now one of these terms m its general form. Let ni,n2, 713,124 G {—2,-1,1,2}, 
ri,r2 e {pi{xi),p2{xi), -p3{xi), ~pi{xi)} and 73,74 G {pi (6)7^2(6), -^3(6): -^4(6)} (with 
TT-i < ??2 < n3 < Hi and ri < r2 < 73 < 74), then wc want to compute 



de2 



G(-,a;2 - i2)d£,2dx2 



e=0 



We can spht up the integral over [ni(5;, +eri, n25h + S72] x [jisiJ;, + £r3, ^4(56 + er4] into nine integrals 
over the domains 



[n25b, n2Sb + e72] x [n3Sb + era, na^fc], 
[712(56, n25b + £72] X [n4(56, n^Sb + er], 
[rtifJb, n2(5fc] X [713(56,714(56], 
[?7.i(56 + eri, 771(56] x [713^6 + era, 773(56], 
[7ii(5b + eri, 771(56] x [7745{,, 774(56 + er4]. 



[772(56, "2(56 + £^2] X [773(56, 714(5(,], 
[71l(5b, 772(56] X [773(5fc + er3,773(5b], 
[77l(5(,, 772(56] X [774(56,774(56 + £7-4], 

[771(56 + eri, 771(56] x [773^6, ?74(56], 



We compute two of these integrals. The others are computed in a similar vein. G,2 denotes the 
partial derivative of G with respect to its second argument. 



de^ 

d 
de 



n2&b+sr2 pnsSb 



r2 rO 



nsSb+era 



G{-,X2 -£.2)d£,2 dX2 



e=0 



— / / e G{-,£{x2~^2) + {ri2-n3)5b)d^2dx2 

no _ 
2£G{-,£{X2 - 6) + (772 - 773)^6) + £^{X2 - 6)G,2(-, £(i2 - 6) + ("2 - 773)^6) 
• 3 

rr2 rO 

G{-, (772 - 773)(5b) d£,2 dx2 

I r3 

-2r2r3G(-, (772 - 773)Jf,). 



(i^2 dX2 



E = 



33 



Another kind of integral we encounter is 



de2 



/ / G{-,X2 - £,2) d^2dX2 

J n2Sb Jn^St 



d_ 

de 



InsSi 
r-2 r{n2—n4^)6iy 



e=0 



eG{-,ex2 + C2) dS,2 dx2 



(ri2-n3)i5b 
r2 /•(n2-n4)<56 



£=0 



'I {n2-n3)Si, 

a(n2-n4)<5b 
2i2G,2{-,i2)dl2dX2 
^n2-n'j,)Si 



G{-,ex2 + ^2) + £X2G^2{-,£S:2 + ^2) dS,2 dx2 



£=0 



= fl \G{{n2 - n3)db) - G{{n2 - ni)5b) ^ 
Combining all integrals we find 

I = -2r2r3G'(-, (7^2 - n:i)5i,) + r\ {g{-, (712 - n:^)^}) - G(-, (^2 - n^E^)^ + 2r2riG{-, {n2 - ni)Sb) 

+ rl (g(-, {112 - n^)5b) - G(-, (nj - 11^)5^)) - rj (g(-, {n2 - 114)6^) - G(-, {ni - ni)5b)^ 

+ 2rir3G(-, (ni - 7i3)(5b) + rj(^G(-, (ni - n4)(56) - G(-, {ni - n3)Sb)^ ~ 2rir4G(-, (ni - n4)(5(,). 

Applying this result to (|38p while keeping in mind that G(— xi,-) = G{xi,-) and G(-,— .T2) = 
G(-, X2), we find 



/E(a;i,6) 



-8p?(a;i) + 8pi(xi)pi(ei) - 2p2(a;i) + 2p2(.Ti)p2(6) 

-8p2(a;i) + 8p3(a^i)P3(Ci) " 2p^(-^i) + 2^4(2:1)^4(^1)] G(xi - 6,0) 

4pi(2;i) - 8pi(a;i)p2(6) + 4^3(2:1) + 4^3(2:1) - 8^3(2:1)^4(6) + 4^4(2:1) G(a;i - 6, (^b) 

+ 8p2(xi) + 16pi(a-i)p3(a) + 8p1(2;i)] G(x-i - ^u2Sb) 

iplixi) + 4pl{xi) + 8p2 (.2:1)^3(6) + ipU^i) + 8pi{xi)p4{£.i) + ipliii)] G{xi - 6, 3<5b) 

2p2(a;i) +4^2(2:1)^4(6) + 2p2(a;i)]G(a;i -6, 4<56), 

where we have used that in ([37]) the integrations over xi and 6 ^-rc indistinguishable. 
Note that, by Corollary [521 for ^ e T^, r e M 



/ G(a; - 6 r-) ^^2; = / G{x,r) dx = -^Ir]. 
Jo Jo 2 

G{q, r) for r G R, we find 



de2 



Using this, as well as equations ([4l[5|) and the equality G{q, r 

'^l\\u, ~ veWl-.^s.^ - -i5bY,{\Mq)\' + \Mq)\') 

+ L^Y1 + 2|P2(<Z)P + 8|p3(g)P + 2\p4{q)\'}G{q, 0) 

- 8{pi{q)P2{q) +ps.{q)pi{q)]G{q, 5b) 

+ l6M<l)Mq)Giq,2Sb) 

- 8{p2{q)p3{q) +Pi{q)P4{q)}G{q,3Sb) 

+4p2{q)Fiiq)G{qASb)\ . (39) 
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Adding the results and ([5^ . we get 



1=0 Jo 



r 2 2 

duv \p'i 



/ 2 , ,2 



^ [{8|pi(g)|2 + 2|p2(g)|2 + S\pM? + 2|p4(g)P}G(g, 0) 

9GZ 



- 8{pi((7)p2(g) +P3(g)p4(9)}G(q,(5&) 

-16pi(g)^(g)G((z,256) 

-8{p2(g)^(g)+Pi(g)^(g)}G(g,35b) 

-4j52(g)^(g)G(g,4(56) 



^5,Y^{\p,{q)\' + \pM? 



Because we have, for all q,q G N, 



27rxq\ , ( 1'Kxq\ 
sm 







L _ 
2 '•^ 



sm 



/ 1'Kxq\ ( 2T:xq\ 
cos 



da; = — 
L 

dx = 0, 



/ 2-Kxq\ ( 2'Kxq\ 
' ' cos ' ' 



-) [-ir) 



dx = Snq, 



the integral over the derivatives in the second variation gives us 



+d„o |(a2,i)^ + (04 J )^ + (^2j)^ + (&4j)^| 



Because p is M'*-valued, p{q) = p{—q)- Furthermore G{—q,X2) = G{q,X2) by equation ([32 
This enables us to write terms as follows, for fc, / G {1, 2, 3, 4}: 

00 

J2Pkiq)Mq)Gi<l.x2) =Pk{0)piiO)GiO,X2) + 2RcJ2M<l)MQ)Giq,x2). 

gGZ q=l 

Note that, for k e {1,2,3,4}, q e Z \ {0}, we have pk{q) = ^ (af' - i^f ^) and pfc(g) = 
^ (^af ^ + ife^'"') and thus, for I £ {1, 2, 3, 4}, 

From the Fourier series (|33p we get furthermore that 

G(0,X2) 



2Vl 



X2\ 



G{q,X2) - ^e-2-l-2|9/L^ for q>l. 
inq 

Using these results in the expression for the second variation yields the desired result. 



D Detailed calculations in the proof of Lemma 14.51 

In this appendix we prove ((22)) . Since < u < 1 we have 3(— 1 + v^) — 4(— 1 + Q log^ f < and 
thus G_ < <^=4> h- > 0. Because |log®i;, the coefficient in front of ('^ in is positive, we 
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know that is positive for C G [0, Ci(^)) U {C2(v), 1], where (1^2 are the u-dcpendent zeros of 
with ^1 < C2- These zeroes are given by 

CiA^) = (81og^w)"^(9- 12u2 + 3i;^ + (41ogi;)(3 + log2u) 

±{225 - 504^2 + 342i;'* - 72v^ + 9v^ + (360 - 288u^ - 72v^) logv 
+ 144 log^ V + (-120 + 96u^ + 24:V^) log^ w - 961og^ v + 161og^ u} ■ 

We take the plus sign in (i and the minus sign in C2. In this way the negativity of (8 log'^ u)^^ 
ensures that Ci ^ C2- Plots of Ci and (2 are given in Figure [TBI 




0.2 0.4 0.6 0.8 

(a) Plot of Ci 




0.2 0.4 0.6 0.8 

(b) Plot of C2 



Figure 16: 

We start by proving that 9 - 12v'^ + + (4 log i') (3 + log^ u) < on (0, 1). The equation 

V— (2v^ - 2v^ + logv) = 8u'' - 4^2 + 1 = 
dv 

has no real solutions on (0, 1) and so ^v-^ (w^ — 2v^ + log^ u + l) = 2u'* — 2u^ + logv < on (0, 1], 
with equality iff u = 1. This in turn shows that 12v-£ (9 - 12^^ + 3v'^ + (41ogu)(3 + log^ v)) = 
w'*-2w2 + log2t; + l > on (0, 1], with equality iff i; = 1, from which 9- 12w2 + 3u'^ + (41ogw)(3 + 
log^ v) < follows. Consequently, since (81og'^u)~^ < 0, we have that ^2 > 0. 
Next we calculate 



— 3 log 1^ - 3 + 
dv 



u^ + 1 



- log^ v] =-- 



12v 



V (u^ + 1)2 



3, 2 
— log V. 

V 



2 (r v^)'^ 

This is equal to zero if and only if log v = which leads to u = e 1+^ 



We will now prove 



v e [0,1] /\v = e 1+-^ 



v = 1. 



Since v ~ I clearly satisfies the equation on the left, it remains to show that there are not more 
solutions. We start by computing 



d 
dv 

dv^ 



e i+" 



e i+" 



On [0, 1] we have 



4u 



[i + v-^Y 

1 

{l + v-^f 



e 1+.. 



1; 



-12z;'' + 8w^ +4) e^+^. 



dv"^ 



e i+- 



-1217"^ + 81)2 + 4 = 



w = 1, 
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showing that ^ I e — i; I has at most one zero on [0, 1] and thus its only zero is at u = 1. 

which in turn shows that also e 1+"^ — v has at most one zero on [0, 1], which is what we set out 
to prove. This now leads us to conclude 

|-(31og.-3+-A_-log3.)^0^. = l. 

This means that 3 log v — 3 + -^jrpi ^ log^ ^ has a minimum at v = 1 and thus this expression is 
positive on (0, 1). Then 

|(9- 12u2 + 3u4 4. (4iogt;)(3 + log2u)) -Slog^wj^ 

- 1 225 - 504^2 + 342u** - 72v^ + 9v^ + (360 - 288^^ - 72v*) \ogv 

+ 144 log^ V + (-120 + 96u^ + 24u'^) log^ v - 96 \og^ v + 16 log^ u| 

= -144(^2 - 1)2 + 144(^4 - 1) logi; - 48(u^ - 1) log^ v 

2 



= 48(u*- 1) [^31ogi;-3 |^1 - ^ ^ j - log^ u 

< 0. (40) 

Note that this also proves that the expression in the square root in C1.2 is positive. Together with 
81og^ V <0 these inequalities give us C2{v) > 1- These results lead to the conclusion that 

G_(u,C) <O^Ce [0,Ci(t^))- 

The other sign possibilities for G_ follow immediately. 
Remark D.l. For the excluded endpoints and 1 we find 

. „ . 5 1 /69 
hmCi = 0, limCi = \ — , 

vio^^ ' ,;Ti^ 2 2 V 5 ' 

,• . 1 , 5 1 /69 

hmC2 = li hniC2 = — I — \ — • 

t^io^ ' vn 2 2V 5 

The limits for v | 1 were found by calculating the first terms in the Taylor expansion of Ci^2- 
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